Observation of breathers in Josephson ladders 
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We report on the observation of spatially-localized excitations in a ladder of small Josephson 
junctions. The excitations are whirling states which persist under a spatially- homogeneous force 
due to the bias current. These states of the ladder are visualized using a low temperature scanning 
laser microscopy. We also compute breather solutions with high accuracy in corresponding model 
equations. The stability analysis of these solutions is used to interpret the measured patterns in the 
I — V characteristics. 
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The present decade has been marked by an intense the- 
oretical research on dynamical localization phenomena in 
spatially discrete systems, namely on discrete breathers 
(DB). These exact solutions of the underlying equations 
of motion are characterized by periodicity in time and 
localization in space. Away from the DB center the sys- 
tem approaches a stable (typically static) equilibrium. 
(For reviews see [0, §). These solutions are robust to 
changes of the equations of motion, exist in translation- 
ally invariant systems and any lattice dimension. Dis- 
crete breathers have been discussed in connection with 
a variety of physical systems such as large molecules, 
molecular crystals ||], spin lattices to name just 

a few. 

For a localized excitation such as a DB, the excita- 
tion of plane waves which might carry the energy away 
from the DB does not occur due to the spatial discrete- 
ness of the system. The discreteness provides a cutoff for 
the wave length of plane waves and thus allows to avoid 
resonances of all temporal DB harmonics with the plane 
waves. The nonlinear ity of the equations of motion is 
needed to allow for the tuning of the DB frequency [^]. 

Though the DB concept was initially developed for 
conservative systems, it can be easily extended to dissipa- 
tive systems [g. There discrete breathers become time- 
periodic spatially localized attractors, competing with 
other (perhaps nonlocal) attractors in phase space. The 
characteristic property of DBs in dissipative systems is 
that these localized excitations are predicted to persist 
under the influence of a spatially homogeneous driving 
force. This is due to the fact, that the driving force com- 
pensates the dissipative losses of the DB. 

So far research in this field was predominantly theo- 
retical. Identifying and analyzing of experimental sys- 
tems for the direct observation of DBs thus becomes a 
very actual and challenging problem. Experiments on 
localization of light propagating in weakly coupled opti- 
cal waveguides |7| , low-dimensional crystals and anti- 
ferromagnetic materials have been recently reported. 




FIG. 1. Schematic view of an annular ladder. Josephson 
junctions are indicated by crosses (x). 

In this work we realize the theoretical proposal to 
observe DB-like localized excitations in arrays of coupled 
Josephson junctions. A Josephson junction is formed be- 
tween two superconducting islands. Each island is char- 
acterized by a macroscopic wave function ^ ~ e^^ of the 
superconducting state. The dynamics of the junction is 
described by the time evolution of the gauge- invariant 
phase difference (f = O2 — Oi — ^ J A ■ ds between adja- 
cent islands. Here $0 is the magnetic flux quantum and 
A is the vector potential of the external magnetic field 
(integration goes from one island to the other one). In the 
following we consider zero magnetic fields A = 0. The 
mechanical analogue of a biased Josephson junction is a 
damped pendulum driven by a constant torque. There 
are two general states in this system: the first state cor- 
responds to a stable equilibrium, and the second one cor- 
responds to a whirling pendulum state. When treated for 
a chain of coupled pendula, the DB corresponds to the 
whirling state of a few adjacent pendula with all other 
pendula performing oscillations around their stable equi- 
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libria. In an array of Josephson junctions the nature of 
the coupling between the junctions is inductive. A local- 
ized excitation in such a system corresponds to a state 
where one (or several) junctions are in the whirling (re- 
sistive) state, with all other junctions performing small 
forced oscillations around their stable equilibria. Accord- 
ing to theoretical predictions [|ll[ , the amplitude of these 
oscillations should decrease exponentially with increasing 
distance from the center of the excitation. 

We have conducted experiments with ladders consist- 
ing of Nb/Al-AlO^c/Nb underdamped Josephson tunnel 
junctions [|l2|. We investigated annular ladders (closed in 
a ring) as well as straight ladders with open boundaries. 
The sketch of an annular ladder is given in Fig. |^. Each 
cell contains 4 small Josephson junctions. The size of 
the hole between the superconducting electrodes which 
form the cell is about 3x3 /im^ . Here we define vertical 
junctions (JJy) as those in the direction of the exter- 
nal bias current, and iiorizontai junctions (JJh) as those 
transverse to the bias. Because of fabrication reasons we 
made the superconducting electrodes quite broad so that 
the distance between two neighboring vertical junctions 
is 30 /im as can be seen in Fig. ^. The ladder voltage 
is read across the vertical junctions. According to the 
Josephson relation, a junction in a whirling state gen- 
erates a dc voltage V = (^^^, where (...) means 
the time average. In order to force junctions into the 
whirling state we used two different types of bias. The 
current Ib was uniformly injected at every node via thin- 
film resistors. Another current /locai was applied locally 
across just one vertical junction. 

We studied ladders with different strength of horizon- 
tal and vertical Josephson coupling determined by the 
junction areas. The ratio of the junction areas is called 
the anisotropy factor and is expressed in terms of the 
junction critical currents r] = /ch/^cV- If this factor is 
equal to zero, vertical junctions will be decoupled and 
can operate independently one from another. Measure- 
ments have been performed at 4.2 K. The number of cells 
N in different ladders varied from 10 to 30. The discrete- 
ness of the ladder is expressed in terms of the parame- 
ter /^L = 27rL/cv/*^0 5 where L is the self-inductance of 
the elementary cell of the ladder. The damping coeffi- 
cient a = ^o/{27rIcCR^) is the same for all junctions 
as their capacitance C and resistance scale with the 
area and Cu/Cy = i^Nv/^NH = V- The damping a in 
the experiment can be controlled by temperature and its 
typical values are between 0.1 and 0.02 . 

We have measured the dc voltage across various verti- 
cal junctions as a function of the currents /locai and Ib- In 
order to generate a localized rotating state in a ladder we 
started with applying the local current /locai > 2/cH+^cV- 
This switches one vertical and the nearest horizontal 
junctions into the resistive state. After that /locai was 
reduced and, simultaneously, the homogeneous bias Ib 



was tuned up. In the final state we kept the bias Ib con- 
stant and reduced /locai to zero. Under these conditions, 
with a spatially-homogeneous bias injection, we observed 
a spatially-localized rotating state with non-zero dc volt- 
age drops on just one or a few vertical junctions. 

Various measured states of the annular ladder in the 
current-voltage Ib — V plane with /locai = are pre- 
sented in Fig. ^ The voltage V is recorded locally on 
the same vertical junction which was initially excited by 
the local current injection. The vertical line on the left 
side corresponds to the superconducting (static) state. 
The rightmost (also the bottom) curve accounts for the 
spatially- homogeneous whirling state (all vertical junc- 
tions rotate synchronously). Its nonlinear Ib{V) shape 
is caused by a strong increase of the normal tunneling 
current at a voltage of about 2.5 mV corresponding to the 
superconducting energy gap. The series of branches rep- 
resent various localized states. These states differ from 
each other by the number of rotating vertical junctions. 
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FIG. 2. Current- volt age Ib—V plane for an annular ladder 
with the parameters A/" = 30, ?7 = 0.44 and /3l = 2.7. Digits 
indicate the number of rotating vertical junctions. 

In order to visualize various rotating states in our lad- 
ders we used the method of low temperature scanning 
It is based on the mapping of a 



laser microscopy |13 



sample voltage response as a function of the position of 
a focused low-power laser beam on its surface. The laser 
beam locally heats the sample and, therefore, introduces 
an additional dissipation in the area of few micrometers 
in diameter. Such a dissipative spot is scanned over the 
sample and the voltage variation at a given bias current 
is recorded as a function of the beam coordinate. The 
resistive junctions of the ladder contribute to the volt- 
age response while the junctions in the superconducting 
state show no response. The power of the laser beam is 
modulated at a frequency of several kHz and the sample 
voltage response is measured using a lock-in technique. 



2 



Several examples of the ladder response are shown in 
Fig. ^as 2D gray scale maps. The spatially- homogeneous 
whirling state is shown in Fig. ||(A). Here all vertical 
junctions are rotating but the horizontal ones are not. 
Fig. H(B) corresponds to the uppermost branch of Fig. |^. 
We observe a localized whirling state expected for a DB, 
namely a rotobreather In this case 2 vertical junc- 
tions and 4 horizontal junctions of the DB are rotating, 
with all others remaining in the superconducting state. 
The same state is shown on an enlarged scale in Fig. |^(C). 
Fig. H(D) illustrates another rotobreather found for the 
next lower branch of Fig. ^ at which 4 vertical junctions 
are in the resistive state. The local current at the begin- 
ning of each experiment is passing through the vertical 
junction being the top one on each map. In Fig. ||(E), 
which accounts for one of the lowest branches, we find 
an even broader localized state. Simultaneously, on the 
opposite side of the ring we observe another DB excited 
spontaneously (without any local current). An interest- 
ing fact is that in experiments with open boundary lad- 
ders (not closed in a ring) we also detected DBs with 
even or odd numbers of whirling vertical junctions. 




FIG. 3. Whirling states measured in the annular ladder 
using the low temperature scanning laser microscope: (A) 
spatially- homogeneous whirling state, (B)-(E) various local- 
ized states corresponding to discrete breathers. 



Various states shown in Fig. g account for different 
branches in the Ib — V plane in Fig. ^ Each resistive 
configuration is found to be stable along its particular 
branch. On a given branch the damping of the junc- 
tions in the rotating state is compensated by the driving 
force of the bias current Ib • The transitions between the 



branches are discontinuous in voltage. In Fig. g, we see 
that all branches of localized states lose their stability 
at a voltage of about 1.4 mV. Furthermore, as indicated 
in the inset on Fig. ^, a peculiar switching occurs: upon 
lowering the bias current Ib the system switches to a 
larger voltage. According to our laser microscope obser- 
vations, the lower is the branch in Fig. |^, the larger is 
the number M of resistive vertical junctions. The slope 
of these branches is dV/dlB ~ MR^y/{M + r^), thus the 
branches become very close to each other for large M. 
The fact that the voltage at the onset of instability is in- 
dependent of the size of the breather, indicates that the 
instability is essentially local in space and occurs at the 
border between the resistive and nonresistive junctions. 

The occurrence of DBs is inherent to our system. We 
have also found various localized states to arise with- 
out any local current. Namely, when biasing the ladder 
by the homogeneous current Ib slightly below NIcv^ we 
sometimes observed the system switching to a spatially- 
inhomogeneous state, similar to that shown in Fig. ||(E). 

To interpret the experimental observations, we ana- 
lyze the equations of motion for our ladders (see Q] 
for details). Denote by (p^,(p^,(f^ the phase differences 
across the Ith vertical junction and its right upper and 
lower horizontal neighbors. Using \/(pi = (pi-\-i — (pi and 
Aipi = (pi-^i + (pi-i — 2(pi, the Josephson equations yield 

+ a^r + siiK^r = 7 - -^(-A^r + v^ti - V(^ti) (i) 

PL 

*.N„^r+si„^f = ^(^|.-^^+v,n (3) 

Here 7 = Ib/{NIcv). First, we compute the dispersion 
relation for Josephson plasmons (f ex e^i^^-^t) at a = 
in the ground state (no resistive junctions). We obtain 
three branches: one degenerated with u = 1 (horizontal 
junctions excited in phase), the second one below 00 = 1 
with weak dispersion (mainly vertical junctions excited) 
and finally the third branch with the strongest dispersion 
above the first two branches (mainly horizontal junctions 
excited out of phase), cf. the inset in Fig. ||. The region of 
experimentally observed breather stability is also shown. 
Note, that for DBs with symmetry between the upper 
and lower horizontal junctions the voltage drop on the 
horizontal junctions is half the drop across the vertical 
ones. This causes the characteristic frequency of the DB 
to be two times smaller than the value expected from the 
measured voltage drop on the vertical junctions O]. 
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FIG. 4. Current- volt age characteristics for numerically 
obtained breather states with 2, 4, 6, and 8 whirling vertical 
junctions {a = 0.07). The inset shows the dispersion relation 
of an annular ladder with 7 = 0.3 and N — 30, 77 = 0.44, 
/3l = 2.7. 

In order to compare experimental results of Fig. ^ 
to the model given by Eqs. (0)-(||) we have integrated 
the latter equations numerically. We also find localized 
breather solutions, in particular solutions similar to the 
ones reported in previous numerical studies jl^. These 
solutions are generated with using initial conditions when 
M vertical junctions of the resistive cluster (cf. Fig. ||) 
have (/9 = and ip — 2V and the horizontal junctions 
adjacent to the vertical resistive cluster have cp = and 
(f = while all other phase space variables are set to 
zero. The obtained current-voltage characteristics are 
shown in Fig. |[ Note that the superconducting gap 
structure and the nonlinearity of slopes are not repro- 
duced in the simulations, as we use a voltage independent 
dissipation constant a in of (|l|)-(^). We find several insta- 
bility windows of DB solutions, separating stable parts 
of the current-voltage characteristics. 

In addition to direct numerical calculation of Ib — 
V curves, we have also computed numerically exact 
breather solutions of (|l|)-(§) by using a generalized New- 
ton map jl|]. We have studied the linear stability of the 
obtained breather [|l[ by solving the associated eigenvalue 
problem. The spatial profile of the eigenmode which 
drives the instability (associated with the edges of the 
instability windows in Fig. 4) is localized on the breather, 
more precisely on the edges of the resistive domain. This 
is in accord with the experimental observation (Fig. 3) 
where several independent breathers can be excited in 
the system. 

The DB states turn to be either invariant under (p^ ^ 
(p^ transformation or not. Both such solutions have been 
obtained numerically. To understand this, we consider 
the equations of motion (||)-(|3|) in the limit r] ^ and 
look for time-periodic localized solutions. In this limit 



the brackets on right hand side of (||) and (||) vanish, and 
vertical junctions decouple from each other. Let us then 
choose one vertical junction with / = to be in a resistive 
state and all the others to be in the superconducting 
state. To satisfy periodicity in the horizontal junction 
dynamics we arrive at the condition 
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and a similar set of choices for (p^i,(f^i, with all other 
horizontal phase differences set to zero. Here k is an ar- 
bitrary positive integer. Continuation to nonzero 77 val- 
ues should be possible [||. The symmetric DBs in Fig. 3 
correspond to k = 2. The mentioned asymmetric DBs 
correspond to /c = 1. The current- voltage characteristics 
for asymmetric k = 1 DBs show a different behavior from 
that discussed above. These DBs are stable down to very 
small current values, and simply disappear upon further 
lowering of the current, so that the system switches from 
a state with finite voltage drop to a pure superconducting 
state with zero voltage drop. 

The observed DB states are clearly different for the well 
known row switching effect in 2D Josephson junction ar- 
rays. The DBs demonstrate localization transverse to the 
bias current (driving force), whereas the switched states 
of non-interacting junction rows are localized along the 
current. At the same time, DB states inherent to Joseph- 
son ladders are closely linked to the recently discovered 
meandering effect in 2D arrays jlj]. 

In summary, we have experimentally detected various 
types of rotobreathers in Josephson ladders and visual- 
ized them with the help of laser microscopy [|l5[. Our 
experiments show that DBs in Josephson ladders may 
occupy several lattice sites and that the number of oc- 
cupied sites may increase at specific instability points. 
The possibility of exciting DBs spontaneously, without 
using any local force, demonstrates their inherent char- 
acter. The observed DBs are stable in a wide frequency 
range. Numerical calculations confirm the reported inter- 
pretation and allow for a detailed study of the observed 
instabilities. 
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